High Rayleigh number turbulent convection in a gas near the gas-liquid critical point. 
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SFg in the vicinity of its critical point was used to study turbulent convection up to exceptionally 
high Rayleigh numbers, Ra, (up to 5 • 10 14 ) and to verify for the first time the generalized scaling 
laws for the heat transport and the large scale circulation velocity as a function of Ra and the 
Prandtl number, Pr, in very wide range of these parameters. The both scaling laws obtained are 
consistent with theoretical predictions by B.Shraiman and E.Siggia, Phys. Rev. A 42, 3650 (1990). 

PACS numbers: 47.27.-i; 47.27.Te; 44.25.+f 



Turbulent convection has recently attracted a lot of at- 
tention due to possibility to tune and control the relevant 
parameters with unbeatable precision. The fundamental 
aspect of the turbulent convection is a competition be- 
tween buoyancy and shear [jlj . An important issue in our 
understanding of the convective turbulence is experimen- 
tal test of theoretical predictions on the Prandtl number, 
Pr, dependence of global characteristics, such as nondi- 
mensional heat transport, Nu(Ra,Pr), and Reynolds 
number of the large scale circulation flow, Re(Ra, Pr). 
The difficulty to measure the Pr dependence in turbulent 
convection arises from the fact that in conventional fluids 
it cannot be varied substantially but only by changing 
the fluid which is not always a simple task. The verifica- 
tion of these scaling relations will be a crucial test of the 
theory. 

In this Letter we present results on the Pr dependence of 
global properties of the flow, namely, global heat trans- 
port, characterized by Nu, and large scale circulation 
velocity, characterized by Re. 

There are two possibilities to reach high Ra number 
convection regime: either to increase temperature dif- 
ference across the cell or to vary physical parameters, 
which appear in the expression for Ra = 9a ^ L , where 
AT is the temperature difference across the cell, g is 
the gravitational constant, a is the thermal expansion, 
and k and v are the thermal diffusivity and kinematic 
viscosity, respectively. The latter possibility was used 
successfully in compressed gases (He and SFq) at almost 
constant value of Pr about one ^-||. Another system 
in which this method of variation of Ra number in a 
wide range can be used, is a gas near its gas-liquid crit- 
ical point(CP). It was realized long time ago that heat 
transport is enhanced dramatically near T c @. But it 
was only recently shown experimentally that the critical 
temperature difference for the convection onset decreases 
drastically with the distance from T c due to strong vari- 
ations of thermodynamical and kinetic properties of a 
gas in the vicinity of CP Q|. Singular behaviour of the 
thermodynamic and kinetic properties of the fluid near 
T c provides the opportunity both to reach extremal val- 



ues of the control parameter Ra and to scan Pr over 
an extremely wide range . All these features make the 
system unique in this respect |]||. 

However, the most exciting aspect of the system is the 
possibility to perform Laser Doppler velocimetry(LDV) 
measurements that we recently demonstrated |^J^] . Small 
temperature differences used to reach high Ra, lead to 
rather small fluctuations of refraction index in the flow. 
It gives us the possibility to use a standard LDV tech- 
nique. On the other hand, the problem of seeding parti- 
cles in a close flow of a gas, which coagulate and sediment 
rather fast, turns out to be unsolvable at the moment. 
Fortunately, we discovered that the existence of the criti- 
cal density fluctuations provided us the possibility to per- 
form LDV measurements of the velocity field in a rather 
wide range of closeness to the CP between 3 ■ 10 -4 and 
1CP 2 in reduced temperature r = (T — T c )/T c , where T 
is the mean cell temperature. The upper limit is defined 
by a small scattering amplitude, and the lower limit- by 
multiscattering in the large convection cell |||9). 
At the same time those obvious advantages come to- 
gether with limitations and new features which should 
be taken into account and studied. Strong dependence 
of gas properties on the closeness to T c manifests itself 
in nonuniform distribution of density in a gravitational 
field and variation of coefficients in the Navier-Stokes 
equation with temperature and density. The former, the 
well-known gravity effect |Tof , can be significant even at 
t ks and 10 cm height cell: the density difference 

across the cell reaches 1% that leads to rather signif- 
icant variations in the fluid properties. Fortunately, 
a temperature gradient compensates the gravity when 
heating from below ||, and strong convection reduces 
the characteristic size of nonuniformity to one of about 
the boundary layer height. This leads to reduction of 
the density nonuniformity to a tolerant level much below 
0.1%. 

In the Boussinesq approximation fluid properties are as- 
sumed to be constant despite the temperature gradient 
across the cell, except for the buoyancy term. Accord- 
ing to [O] the degree of deviation from the Boussinesq 
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approximation is adequately described by the ratio of 
the temperature drop across the top boundary layer to 
the temperature drop across the bottom boundary layer. 
As the experiment shows [ fi"T| , significant deviations in 
the measurements of global transport properties of the 
turbulent convection occurs if this ratio reduces below 
0.5. According to our estimates measurements at Ra up 
to 10 15 can be done while non-Boussincsq effects are still 
relatively small. In the data presented the temperature 
drop ratio was above 0.7. So small deviations from the 
scaling behaviour were observed at highest Pr and Ra. 
There are another aspects, related to the proximity to 
T c such as compressibility (besides adiabatic gradient) 
and breaking down of hydrodynamic description due 
to macroscopical size of thermodynamical fluctuations. 
Simple estimates show that both these factors do not 
play any significant role in the range of parameters un- 
der studies ||. 

Contrary to the compressed gas convection f|-[|], here 
one cannot cover wide range of Ra for one value of Pr. 
On the other hand, the data presented cover the range 
of Ra from 10 10 at low values of Pr (Pr = 8) far from 
T c and AT = YlmK up to the largest attainable in a 
laboratory Ra — 5 • 10 14 at high values of Pr close to 
T c . The whole range of the reduced mean temperature 
was 4 • 10~ 2 > t > 2 ■ 10~ 4 . And finally, as a result of 
large compressibility and relatively large cell height the 
adiabatic temperature gradient was observed |T^ |. This 
effect will be discussed below. 

The experiment we present here, was done with a high 
purity gas SF 6 (99.998%) in the vicinity of T c and at 
the critical density (p c = 730kg /m 3 ). This fluid was 
chosen due to the relatively low critical temperature 
(T c = 318.73 K) and pressure(P c = 37.7 bar) and well- 
known thermodynamic and kinetic properties far away 
and in the vicinity of CP. This gas was widely used to 
study the equilibrium critical phenomena. The cell is a 
box of a cross-section 76x76 mm 2 formed by 4 mm plex- 
iglass walls, which are sandwiched between a Ni-plated 
mirror-polished copper bottom plate and a 19 mm thick 
sapphire top plate of L = 105 mm apart. 
The cell is placed inside the pressure vessel with two side 
thick plastic windows to withstand the pressure differ- 
ence up to 100 bar. So the cell had optical accesses from 
above through the sapphire window and from the sides. 
They were used for both shadowgraph flow visualization 
and for LDV. The pressure vessel was placed inside a wa- 
ter bath which was stabilized with rms of temperature 
fluctuations at the level of 0.4 mK. The gas pressure 
was continuously measured with 1 mbar resolution by 
the absolute pressure gauge. Together with calibrated 
100f2 platinum resistor thermometer they provide us the 
thermodynamic scale to define the critical parameters of 
the fluid and then to use a parametric equation of state 
developed recently for SF 6 |0 . 

We should point out here that it is not obvious at all 



that the critical parameters, defined at equilibrium con- 
ditions, can be used for a strongly non-equilibrium state 
to define closeness to CP. Contrary, as suggested by the- 
ory p3 and experiment p5| , turbulent flow in a binary 
mixture near its consolute CP can suppress the critical 
concentration fluctuations. At Re comparable with that 
achieved in our experiment, Pine et al p5[ observed the 
critical temperature depression up to 50 mK. We used 
the enhancement of the heat transport (or decrease in the 
temperature difference across the cell at the fixed heat 
flux) near CP as an indicator of the closeness to CP. 
The critical pressure and temperature obtained by this 
procedure agree within ±10mK with those defined inde- 
pendently by light scattering at equilibrium conditions 

The heat transport measurements were conducted at 
the fixed heat flux and temperature of the top sapphire 
plate, while the bottom temperature was measured by a 
calibrated thermistor epoxied in it. Local temperature 
measurements in a gas were made by three 125/im ther- 
mistors suspended on glass fibers in the interior of the 
cell (one at the center, and two about half way from the 
wall). Local vertical component velocity measurements 
at about L/4 from the bottom plate were conducted 
by using LDV on the critical density fluctuations 
Shadowgraph visualization was used mostly to get qual- 
itative information about structures and characteristic 
time and length scales in the flow mostly of the top and 
bottom boundary layers. 

From the measurements of the heat transport and, par- 
ticularly, of the velocity we realized that at a small but 
finite temperature differences, much larger than defined 
by the critical Ra for convection onset, there exists a 
mechanically stable state. This temperature difference 
AT a d is defined by the adiabatic temperature gradient for 
the convection onset. The latter results from a fluid com- 
pressibility and is defined as AT aa -/L — gTaCp^ 1 [Q, 
where C p is the heat capacity at constant pressure. As 
follows from simple thermodynamics and critical diver- 
gences of thermodynamical parameters as T approaches 
T c , the adiabatic temperature gradient saturates at a 
finite value AT ad /L = gp(dP/dT)- 1 , which for SF 6 and 
L = 10.5 cm gives AT ad = 9.5mK for r < 10~ 2 , So with 
available temperature stability and resolution it can be 
measured. Then, as shown in p), Ra should be modified 

The most sensitive probe to detect the convection onset 
in our case was the local velocity measurements. The re- 
sults of the mean vertical velocity V m and rms of vertical 
velocity fluctuations V s measurements at r = 8 ■ 10~ 4 are 
presented in Fig.l. Both signals show clear transition at 
AT a d — 9.5 ± 0.5mK , that agrees well with the theo- 
retical value |||l(J. These measurements were done for 
several values of the reduced temperature in the range 
8 ■ 10~ 3 > t > 2 ■ 10~ 4 . It was found that AT ad (r) is con- 
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stant and independent of r in agreement with the theory 
M. This value of AT a( j was used further to correct the 
Ra values. As a result of this correction for each Pr both 
Nu vs Aa(Fig.2) and V s vs Aa(Fig.3) dependences show 
scaling laws in rather wide range of Ra 
The heat transport measurements were done in the range 
1.6 • 10~ 2 > t > 2 • 10~ 4 . In order to keep Pr constant, 
t was kept constant, while AT changed during the mea- 
surements. The range of r was limited from higher 
temperature end by mechanical stability of construct- 
ing materials (mostly plexiglass) and from lower end 
by the temperature stabilization of the system. Heat 
transport measurements were also conducted far away 
from CP at P = 20bar, T = 303A and the density 
p = 0.18(;/cm 3 , that corresponds to Pr = 0.9, and at 
P = 50bar, T = 323A and p = 1.07g/cm 3 , that cor- 
responds to Pr = 1.5. The data far away from CP 
cover the range of Ra between 10 9 and 5 • 10 12 with 
the temperature differences across the cell from about 
0.1 A" till 10K. By making all appropriate corrections 
for heat losses through the lateral walls, gas and insula- 
tion outside the cell, and for a temperature drop across 
the top sapphire plate, we found that the data on the 
heat transport far away from CP are in a good agree- 
ment with a 2/7 law for 10 9 < Ra < 5 • 10 12 [§. These 
data combined with the heat transport measurements 
in the same cell in the vicinity of CP fot higher val- 
ues of Pr can be scaled and presented in a power law 
form: Nu = o.22Aa - 3±0 03 Pr-°- 2±004 (Fig.4). Just a 
few points at the highest Pr and Ra deviate from this 
scaling due to the non-Boussincsq effect. This scaling is 
consistent with predictions of ref jl| . 
We visualized a structure of the large scale flow by us- 
ing shadowgraph technique through the top and side 
windows. Large diameter with narrow focal width imag- 
ing optics enabled us to visualize narrow slices across the 
cell and scan it in both vertical and horizontal directions. 
The flow was recorded at various Ra and Pr. This vi- 
sualization confirmed the picture of up- and down-going 
circulation jet flow along the cell diagonal, which forms 
the top and bottom turbulent boundary layers. The bulk 
of the cell appears homogeneous. The turbulent charac- 
ter of the flow at the boundary layers is seen as a rapid 
and "violent" horizontal motions ||||. 
The circulation frequency which corresponds to the 
travel time of a fluid element passing through one cycle 
of the large eddy, is observed as a peak in the power 
spectra for both the velocity and temperature fluctua- 
tions (see the insert in Fig.5). We also measured directly 
by LDV the large scale circulation velocity. However, 
the best results were obtained by extracting the peak 
frequency from the velocity power spectra at various 
Ra and Pr which can be presented in a power law 
form: Re = 2.6Ra 0A3±a02 Pr~°- 75±Q02 (Fig.5). Here 
Re = Af p L 2 /v, and f p is the peak frequency. The data 
from the low frequency peak in the temperature spectra 



far away and close to CP and from LDV measurements 
are consistent with this scaling law. We would like to 
point out that Re reaches values up to 10 5 at the highest 
values of Ra. 

The dependence of Re on Pr has almost the same expo- 
nent (within the measurement uncertainty) as one found 
recently in convection in helium for much narrow range 
of Pr |l7j ], and agrees rather well with the theoretical 
prediction 5/7 in ref Q. However, the scaling of Re with 
Ra differs significantly and close to 3/7 rather to 0.5 Jl7| . 
We can also verify consistency of the exponents of the 
both scaling laws obtained by using one of them, e.g. 
for i?e, and the exact relation for the dissipation in a 
bulk turbulent regime PrRa(Nu-l) ~ Re 3 Pr 3 0. The 
resulting scaling relation Nu ~ Ra 29 Pr~ a25 is con- 
sistent with that found experimentally. Together with 
the unified scaling law in the whole range of Ra and Pr 
under studies for both transport mechanisms it suggests 
that a single mechanism is responsible for these scaling 
laws. Moreover, our data did not indicate any signature 
of the transition to the asymptotic regime in the heat 
transport, discovered recently at relatively low Pr and 
Ra > 10 11 . There are possible two explanations to this 
fact: either at higher Pr the transition occurs at higher 
Ra, or the scatter in the data for different Pr does not 
allow us to observe the transition. 

In conclusion, we present for the first time the general- 
ized scaling laws of the heat transport and the large scale 
flow in respect to Ra and Pr in a wide range of variation 
of there parameters. Both scaling laws are consistent 
with the theoretical predictions of Q. We believe that 
in the light of the recent theoretical reconsiderations of 
these scaling laws [|l8| these data provide new insight on 
the turbulent convection. 
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FIG. 1. LDV measurements of the mean V m (squares) and 
rms V s { circles) vertical velocity component. 

FIG. 2. Nu vs Ra: open squares-uncorrected, and full 
squares-corrected data at Pr = 27. 

FIG. 3. V s vs Ra: circles-uncorrected, and triangles- cor- 
rected data at Pr — 93. 



FIG. 4. Nu vs Ra for different Pr. 



FIG. 5. Re of the large scale circulation vs Ra for dif- 
ferent Pr: circles- 27; up-triangles- 45; down-triangles- 93, 
diamonds- 190. The insert: power spectra density of vertical 
velocity fluctuations at Ra — 3 ■ 10 12 at Pr — 93. 
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